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ABSTRACT 

Consider  the  model  of  two  competing  species  with  non-linear  diffusion  and  three  types  of  functional 
dependencies.  The  first  type  of  dependence  corresponds  to  Malthusian  demographic  processes,  the  second  -  ferhulst 
(logistic  population),  and  the  third  -  the  population  of  "Allee".  The  common  element  to  this  kind  of  description  is  the 
presence  of  a  linear  source,  and  in  the  descriptions  of  the  populations  of  Ferhulst  and  Allee  type  are  also  nonlinear  sinks. 

Offered  suitable  initial  approximations  for  quickly  convergence  iterative  process.  Numerical  experiments  are 
carried  out  with  visualization  for  different  values  of  parameters  of  the  system  of  biological  population  of 
Kolmogorov-Fisher  type. 

KEYWORDS:  Biological  Population,  Nonlinear  System  of  Partial  Differential  Equations,  Initial  Approximation, 
Numerical,  Iterative  Process,  Self-Similar  Solutions 

INTRODUCTION 

Blow-up  regimes  in  spatial  localization  in  open  dissipative  systems  described  by  models  with  nonlinear  diffusion 
[1-2].  Some  cross-diffusion  models  discussed  in  several  papers  in  physical  systems  (plasma  physics)  [2,3],  in  chemical 
systems  (dynamics  of  electrolyte  solutions  [4]),  biological  systems  (cross-diffusive  transport  [5],  the  dynamics  of 
population  systems  [6,  7-8],  in  ecology  (return  dynamics  of  forest  structure  [9,10]),  in  seismology  model  Burridj-Knopoff 
describing  the  interaction  of  tectonic  plates  [11-12].  In  the  last  ten  years  in  the  study  of  the  growth  and  development  of 
tumors  actively  used  mathematical  models  with  cross-diffusion  [13-14]  (as  well  as  the  system  of 
reaction-diffusion-advection  [14]). 

Let  us  explain  what  we  mean  by  the  term  "cross-diffusion"  (or  cross-diffusion).  Consider  the  following  system  of 
two  equations  in  one-dimensional  case: 


=  f(u1,u2)  +  Dl 


d2u,     ,  d 

 2    +  kl   

dx  dx 


dux 
~dt 

Bv      ,        .     _  d2u2  d 

—  =  g(uvU2) +  02-^-2-+^ — 


du2^ 


Q1(ul,u2) 


V 


dx 
)«, 
dx 


(D 


At  h\=h2=0  mathematical  model  (1)  is  a  system  of  reaction-diffusion  with  the  diffusion  coefficients 
Dj  >  0,  D2  >  0  (at  least  one  Di  ^  0 ).  In  the  case  where  at  least  one  of  the  coefficients  hi  ^  0  (_sign  can  be  any), 
system  (1)  is  a  cross-diffusion. 
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Linear  cross-diffusion  corresponds  to  Qi(ul,U2)  =  const  for  i=l,2;  for  nonlinear 
cross-diffusion-  Qi(ul,U2)  ^  const  at  least  one  i. 

Cross-diffusion  means  a  spatial  displacement  of  an  object  described  by  one  of  the  variables  takes  place  by 
diffusion  of  another  object,  the  other  variable  described.  At  the  population  level,  the  simplest  example  of  a  parasite 
(the  first  item),  located  inside  the  «host»  (the  second  object  moves  through  the  diffusion  of  the  owner). 
The  term  «self-diffusion»  (diffuse,  direct  diffusion,  ordinary  diffusion)  involves  the  movement  of  individuals  at  the 
expense  of  the  diffusion  flow  from  areas  of  high  concentration,  particularly  in  the  area  of  low  concentration.  The  term 
"cross-diffusion"  refers  to  the  movement  /  flow  of  individuals  of  one  species  /  substances  due  to  the  presence  of  the 
gradient  of  other  individuals  /  agents.  The  magnitude  of  cross-diffusion  coefficient  can  be  positive,  negative  or  zero. 
Positive  cross-diffusion  coefficient  indicates  that  movement  of  individuals  is  in  a  direction  other  species  of  low 
concentration  occurs  in  the  direction  of  high  concentrations  of  other  species  of  animals  /  substances.  In  nature,  systems 
with  cross  diffusion  are  quite  common  and  play  an  important  role,  especially  in  biophysical  and  biomedical  systems. 

Equation  (1)  is  a  generalization  of  the  simplest  diffusion  model  for  the  logistic  model  of  population  growth 
[16-20]  in  Malthusian  type  ( fl(u]  ,U2)  =  Ux ,  fx(ux,U2)  =  U2 ,  f2(ul,U2)  =  Ul ,  f2(ul,U2)  =  U2),  in  Ferhulst  type 
(fl(ul,U2)  =  Ul(l-U2),fl(ul,U2)  =  U2(l-Ul),f2(u],u2)^u1(l-u2),f2(ul,U2)  =  U2(l-Ul)),  and  Allee  type 
(fl(ux,u2)  =  ux{\-u2>),f{(ul,u2)  =  u2(\-ulh),  f2(ux,u2)  =  ux(\-u2l) ,  f2(ux,u2)  =  u2(\-uf2), 
/?,  >1,  P2  >1)  for  the  case  of  double  nonlinear  diffusion. 

In  case,  when  /3{>\,  f$2  >  1 ,  it  can  be  viewed  also  as  the  equation  of  nonlinear  filtration,  thermal  conductivity 

while  the  impact  source  and  absorption  capacity  which  are  respectively  Ux  ,—U2 ,  U2,—uf2 

Consider  a  spatial  analogue  of  Volterra-Lotka  competition  system  with  a  nonlinear  power  dependence  of  the 
diffusion  coefficient  of  density  of  population.  In  the  case  of  the  simplest  volterra  competitive  interactions  between 
populations  can  be  constructed  numerically,  and  in  some  cases  analytically  heterogeneous  in  space  solutions  [21]. 

STATEMENT  OF  THE  TASK 

Consider  in  the  domain  Q={(t,x):  0<  t  <  °°,  xeR2}  parabolic  system  of  two  quasilinear  equations  of 
reaction-diffusion  task  of  biological  populations  of  Kolmogorov-Fisher  type. 


that  describes  the  process  of  biological  populations  in  a  nonlinear  two-component  environment,  the  diffusion 


(2) 


2  f=0 


—  u20(x) , 
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coefficient  is   equal   to    D^U^1    and    D^u"2 ,    <JX ,  <72 ,  f3x ,  J32    -  positive  real   numbers,    Ux  =  >  0  , 

U2  =  U2(t,x)  >  0  -  required  solutions. 

The  Cauchy  problem  and  boundary  problems  for  the  system  (1)  in  the  univariate  and  multivariate  cases 
investigated  by  many  authors  [17-23]. 

The  aim  of  this  work  is  to  study  the  qualitative  properties  of  solutions  of  the  problem  (2)  on  the  basis  of  the 
self-analysis  and  numerical  solutions  using  the  methods  of  modern  computer  technology,  research  methods  of  linearization 
to  the  convergence  of  the  iterative  process  with  further  visualization.  Found  estimates  for  solutions  and  the  resulting  free 
boundary,  that  gives  the  chance  to  choose  the  appropriate  initial  approximation  [17]  for  each  of  the  values  of  numerical 
parameters. 

Let's  build  system  of  equations  (l)-(2)  -  more  simple  system  of  equations  for  research. 
CONSTRUCTION  SELF-SIMILAR  SYSTEM  OF  EQUATIONS 

Let's  construct  a  system  of  equations  by  the  method  of  nonlinear  splitting  [3]. 
Replacing  in  (2) 

ux{t,x)  =  ekl\(t,x), 
u2(t,x)  =  ek2'v2(t,x) 

lead  (2)  to  the  form: 


dVj  _  d 

dr  dx 
dv2  _  d 
dl  dx 


a,  OV2 


dx 
ox 


-k^e 


-k2(t)e 


(Ah-°2h-h+h)t,,A,, 


(3) 


Vi|«=o=Vlo(*).  V2|f=0=V20(*)> 

By  choosing  G^k^  ~*~  ^2  ~~  ^1  =  ^2^1  +  kl  —k2,  we  obtain  the  following  system  of  equations, 


__  d 

~dr 

dx 

< 

dv2 

_  d 

dT 

dx 

dx  j 
dv  ^ 


-a2(t)zb^-v2 


dx 


(4) 


where  ax  =  fc,  {a,k2  +k2~k,  f  ,      bx  =  ^2     <*A     ^2+K  ^ 


<7xk2  +  k2  kx 
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t>2  = 


If      =  0  ,  i  —  1,2  ,  the  system  has  the  form: 


3v,  _  3  f 


dr  dx 

3v, 


v 


-fli  (t)vxv2 


dt  dx 


OX 


V 


The  Cauchy  problem  for  system  (4)  in  case,  when  bl  =  b2  =  0  studied  in  the  works  [18,  21]  and  proved  the 
existence  of  a  wave  of  global  solutions  and  blow-up  solutions. 

Below  we  will  describe  one  way  of  getting  self-similar  system  for  the  system  of  equations  (4).  It  consists  in  the 
following.  Find  the  first  solution  of  a  system  of  ordinary  differential  equations 


dv 


— L  =  -a1v1v2A 
dr 

dv2  a 
.  dr 

in  the  form 

v^T)  =  Cl(T  +  TQyK  ,  v2(T)  =  c2fr+Toy* ,  T0  >  0, 

where 


1  1  1 

Ci=l,7i=-^-.c2=l,  y2=  — 
Pi  A 


And  then  the  solution  of  system  (3)-(4)  is  sought  in  the  form 


V1(t,x)=Vi(t)wl(T,x), 

v2 (t, x)  =  vi(t)w2 (t, x)  , 

and  X  =  X(?)  selected  as 


T1(lT)=jv1Cri(^ 


1 


ln(r  +  r), 


(5) 


Then  for  W,  (T,  x),    i  =  1,2  we  get  following  system  of  equation 
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dT 

dx 

dw2 

__  a 

dx 

6  a  d\V7  a 

(Aw?— ^)  +  ^1(w1w2A-w1) 
dx 

OX 


(6) 


where 


(1-^)7 


i/  l-Yx^x  >0, 


(7) 


sf  i-w=o. 

Representation  of  system  (1)  in  the  form  (5)  suggests  that,  when 

« aw0. 


a 
= — 

< 

r  — >  oo ,  ^ ,  — >  0  and 

ar 

dw2 

dx 

_  a 

[dT 

dx 

dx 

a  aw, 

dx 


(8) 


Therefore,  the  solution  of  system  (2)  with  the  conditions  (5)  tends  to  the  solution  of  the  system  (8). 

Let  JxOx  >  1 ,  Y\°\  =  Yiai .  °2  (^2  +  1)  +  A  (^i  +  1)  =  f  i  (^i  +  1)  +  A  (fc2  +  1)  '  C,  >  0  .  In  this  case, 
by  assuming  in  (6) 

wMtU)=yM  <Hxi/<2,  i=i,2, 

and  considering,  that  the  equation  W;(T,x)  without  the  younger  members  always  has  a  self-similar  solution  in 
the  case  of  1  —  y1<71  ^  0,  we  obtain  a  system 

dg       dg     2  dg 

d  (y?%)4%+/^2-y2  yf2)  =  o, 


J£    2   J£     2  dZ 


where  Mi 


i  -  m 


The  study  of  qualitative  properties  of  the  system  (2)  is  allowed  to  perform  the  numerical  experiment  depending  on 
the  values  included  in  the  system  of  numerical  parameters.  For  this  purpose,  as  the  initial  approximation  was  used 
constructed  asymptotic  solutions. 
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The  numerical  solution  of  the  problem  for  linearization  of  system  (2)  were  used  linearization  methods  of  Newton 
and  Picard.  To  build  self-similar  system  of  equations  of  biological  population  used  the  method  of  nonlinear  splitting 
[18,21]. 

NUMERICAL  EXPERIMENT 

To  solve  the  task  (l)-(2)  by  numerical  method  we  construct  a  uniform  grid 

G)h=\xi=ih,  h>0,  i  =  0,l,...,n,  hn  =  l\ 
and  temporary  grid 

,  =  kj  =  Jh\  >  K  >  °>  j  =  QX-,  n,  vn  =  T\. 

Replace  the  task  (l)-(2)  with  implicit  difference  scheme  and  receive  differential  task  with  the  error  o(h2  +  hl). 

As  it  is  known,  the  main  problem  for  the  numerical  solution  of  nonlinear  problems  is  the  suitable  choice  of  the 
initial  approximation  and  the  method  of  linearization  of  system  (2). 

Let's  consider  the  function: 

vl0(t,x)=v1(t)-{a-£2f, 
v20(t,x)=v2(t)-(a-g2)l2, 

where  Vx(t)  =  ek,Vl(t)  and  V2(t)  =  ektV2(t)  functions  defined  above, 

(2  +  a,) 


fa+ifot  +i)-i' 

(2  +  cr2) 
(<rI+lX<72+l)-l 


Record  (a)+  means  =  max(0,a)  .  These  functions  have  the  property  of  finite  propagation  speed  [18,  21]. 

Therefore,  in  the  numerical  solution  of  task  (2)  when   f3x  >  <JX    as  initial  approximation  proposed  functions 

vi0(t,x)  ,i  =  l,2. 

Created  on  input  language  Math  Cad  program  allows  you  to  visually  trace  the  evolution  process  for  different 
values  of  the  parameters  and  data. 
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Numerical  calculations  show  that  in  the  case  of  arbitrary  values  £7  >  0,  (5  >  0  qualitative  properties  of  solutions 
do  not  change.  Listed  below  are  the  results  of  numerical  experiments  for  different  values  of  parameters  (Figure  1  -  3). 


Values  of 
Parameters 


£7  =  1. 1 

kx  =  0.1 
k2  =  0.25 
eps  =  10"3 


£7  =  1.3 

=5 

k2=l 

eps  =  10" 


Results  of  Computational  Experiment 


20  *B 


Figure  1:  Dynamics  of  Malthusian  Population 


Values  of 
Parameters 


£7  =  1.1 

A=i,  K=i 

J32=l,  k2=9 
eps  =  10"3 


£7  =  1.3 

A=l,  K=3 
J32=l,  k2=2 
eps  =  10"3 


Results  of  Computational  Experiment 


Figure  2:  Dynamics  of  the  Logistic  Population 
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Values  of 
Parameters 

Results  of  Computational  Experiment 

(7  =  1.1 

A  =7,  kx=l 
f$2=9,  k2  =5 
eps  =  10"3 

0  5- 

m 

0                  10                 20  30 

C7  =  1.3 

A  =3,  k,=5 
fi2  =  7,  k2  =  3 
eps  =  10"3 

0  5- 

-oS*£ 

o            io           la  ^0 

Figure  3:  Dynamics  of  Allee  Effect 


CONCLUSIONS 

Thus,  it  can  be  expected  that  future  theoretical  and  experimental  research  of  excitable  systems  with 
cross-diffusion  will  make  a  significant  contribution  to  the  study  of  phenomena  of  self-organization  in  all  of  nonlinear 
systems  from  the  micro  and  astrophysical  systems  to  public  social  systems. 
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